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In this paper we study the effect of, well-known, higher order corrections to the allowed beta 
decay spectrum on the determination of anti-neutrino spectra resulting from the decays of fission 
fragments. In particular, we try to estimate the associated theory errors and find that induced 
currents like weak magnetism may ultimately limit our ability to improve the current accuracy and 
under certain circumstance could even largely increase the theoretical errors. We also perform a 
critical evaluation of the errors associated with our method to extract the anti-neutrino spectrum 
using synthetic beta spectra. It turns out, that a fit using only virtual beta branches with a judicious 
choice of the effective nuclear charge provides results with a minimal bias. We apply this method 
to actual data for 235 U, 239 Pu and 241 Pu and confirm, within errors, recent results, which indicate 
a net 3% upward shift in energy averaged anti-neutrino fluxes. However, we also find significant 
shape differences which can, in principle, be tested by high statistics anti-neutrino data samples. 
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I. INTRODUCTION 



Neutrino physics as a branch of experimental physics started with Cowan's and Reines' detection of neutrinos 
produced in a nuclear reactor [1]. Since then, reactor neutrino experiments have played a crucial role in shaping 
our understanding of the physical properties of the neutrino^] for a review see e.g. Ref. [2], and will continue to do 
so in the future [3H5]. In a nuclear reactor there are about 6 /3-decays, and hence neutrinos, per fission or about 
2 x 10 20 neutrinos per second per GW of thermal power. Fortunately, there are only four isotopes whose fission 
make up more than 99% of all reactor neutrinos with an energy above the inverse /3-decay threshold: 235 U, 239 Pu, 
241 Pu and 238 U. Nonetheless, the resulting neutrino flux is a superposition of thousands of /3-decay branches of the 
fission fragments of those four isotopes and thus, a first principles calculation is challenging, even with modern nuclear 
structure data 6 . Therefore, reactor neutrino fluxes from the thermal fission of 235 U, 239 Pu and 241 Pu are obtained by 
inverting measured total /3-spectra which have been obtained in the 1980s at the Institut Laue-Langevin (ILL) [7HS]- 
There are, currently, no corresponding total /3-spectra for 238 U since this isotope is fissioned only by fast neutrons. 
Modern multi-baseline experiments [SHS] are designed to be independent of a very precise flux knowledge at the source, 
however many of the previous measurements employed only one baseline [2] and initial data from Double Chooz [3] 
and Daya Bay g] most likely will be single baseline data, too. 

Recently, a reevaluation of reactor neutrino spectra has been performed [B] and a significant upward shift in the 
predicted fluxes of about 3% was found. This in turn translates into a weakening of existing reactor bounds on the 
absence of 9 e disappearance with wide ranging consequences on the possible existence of a sterile neutrino with a 
Am 2 > 1 eV 2 [ini E] . Here, we present an independent inversion of the ILL /3-spectra into neutrino spectra for the 
isotopes 235 U, 239 Pu and 241 Pu. To this end, we first review various corrections to the allowed /3-decay shape and 
estimate the associated theory errors in section [TTJ We will use virtual /3-branches and apply the corrections obtained 
in section |n] individually for each branch. We then use synthetic data sets and Monte Carlo simulations to quantify 
the bias and statistical errors associated with the inversion procedure in section |III A| In section |III B| we compute 
the effective nuclear charge Z and provide error estimates. We then apply the method which we have developed to 
the actual /3-data and present our result on the neutrino fluxes in section IV In the following section [V] we perform 



a critical comparison of our result with previously obtained neutrino fluxes and point out a possible explanation for 
the flux shift in terms of known nuclear physics. Finally, we conclude in section |VT] 

II. BETA SPECTRA 

In this section we describe corrections to the allowed /3-spectrum shape and provide an estimate of the associated 
uncertainty. Here, we mainly follow the notation of Ref. |12) . The allowed /3-spectrum is given by [13] 

Np(W) = Kp 2 (W-W Q f F(Z,W), (1) 

^ v X 

phase space 

where W — E/(m e c 2 ) + 1 and Wq is the value of W at the endpoint. K is a normalization constant. F(Z, W) is the 
so called Fermi function and given by 

F(Z, W) = 2( 7 + 1){2p R)*to-%™™/P with 7 = 0^^, (2) 

where R the nuclear radius, which in itself has a dependence on A. The Fermi function accounts for the fact that 
the outgoing electron is moving in the Coulomb field of the nucleus and is derived from the solution of the Dirac 
equation for a point-like and infinitely heavy nucleus. The Fermi function is the leading order QED correction to 
nuclear /3-decay. 

In equation [T] the /3-spectrum shape for allowed decays is given, however fission fragments contain a significant 
fraction of decay branches of unique and non-unique forbidden type. For forbidden decays the simple phase space 
factor p 2 (W - W ) 2 is multiplied by [14] 

p 2 + p 2 1 st unique , 2 nd non-unique 
p* + —p^p 2 + p A 2 nd unique , 3 rd non-unique 
pi + Iplp 2 + Ip^p 2 + p 6 3 rd unique , 4 th non-unique . (3) 



1 Throughout this paper we will refer to both neutrino and anti-neutrino simply as neutrino and rely on the context for disambiguation. 
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We see, that the forbidden space phase factor is symmetric under the exchange of p and p v and thus, if we had 
m„ = m e , the whole expression would be symmetric between neutrino and electron energies, too. If we keep the 
measured /3-spectrum fixed and fit a /3-shape to it, any change to the /3-shape which is symmetric between E e and E v 
would not change the neutrino spectrum. This the reason the corrections from forbidden decays to inverted neutrino 
spectra have been found to be small [3-9: despite the overall large contribution of forbidden decays in fission fragments. 

It turns out, that for precision studies a number of additional effects have to be taken into account and the modified 
/3-spectrum becomes |15j 

Nfs(W) = Kp 2 {W - Wq) 2 F{Z, W) L {Z, W) C{Z, W) S{Z, W) Gp{Z, W) (1 + 5 WM W) . (4) 

The neutrino spectrum is obtained by the replacements W — > Wq — W and Gp G v . In the following we will describe 
the physical origin and the actual expression used for each of these additional terms. The size of the contribution of 
each of these terms for a typical /3-decay is shown in figure [l] 



Finite size corrections 



There are two effects from the finite size of the nucleus: the electric charge distribution is no longer point-like and 
the hypercharge distribution is no longer point-like, i.e. the nucleon is moving inside the nuclear potential. For all 
finite size corrections we need to be able to determine the nuclear radiu^] as a function of A and we use the so called 
Elton formula [T7] 

R = 0.0029A 1/3 + 0.0063A~ 1/3 - 0.017A" 1 (5) 

in units of m e c 2 . 

The electromagnetic finite size effect is expressed by Lq and we use the approximation given in Ref. [15] 



L (Z, W) = 1+13^ -WRaZ -aZR 1 1} +<*-! a n (WRr+0Al(R-0.0W4)(aZ) 4 - 5 



(6) 



where the a„ are given by 

6 

a n = Y,b x {aZ) x , (7) 

x=l 

and the b x are taken from table 1 of Ref. [TS] and reproduced in table [V] in appendix Aj We checked that the error 
of this approximation in comparison with numerically exact results given in |14j is smaller than 10~ 5 . Neglecting the 
last three terms in the above expression for Lq still produces accurate results and this form of Lq produces results 
very close to the ones summarized as finite size effects in equation 8 of Ref. [6] . 

The weak interaction finite size correction is of the following form for Gamow- Teller decayfj^] 15 

C(Z, W) = 1 + Cq + C 1 W + C 2 W 2 with (8) 

630 5 35 

C 1 = -'^RaZ+ Vq-R 2 , 
35 9 

In Ref. [TH] it is estimated that the errors on both these corrections, Lq and C, are smaller than 10% of their size. 
Thus, we can neglect this as a relevant source of error. 



2 Differences due to different shapes of the charge distribution, e.g. constant within the nucleus or all charge on the surface are very 
small and can be neglected, see e.g. 1161 . Also, reasonable variations of the Elton formula have only a minor impact on our result, i.e. 
a variation of 10% in A will change the total flux by 0.05%. 

3 Almost all decays relevant in our problem are of this type and the corresponding expression for Fermi decays is numerically very close. 
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Screening correction 



S(Z, W) accounts for screening of the nuclear charge by all the electrons in the atomic bound state and it effectively 
reduces the charge "seen" be the outgoing electron. Our description of the screening correction S(Z, W) follows 
Ref. [TB] and we introduce the following quantities 

r= aZW aZW 

W = W-V 0l v = VW 2 - 1 , y = ^- y= — - Z = Z-l. (9) 

P P 

Vq is the so called screening potential and derived from numerical calculations and can be parametrized as 

V = a 2 Z 4 / 3 N(Z) , (10) 



VI 



and N(Z) is taken as linear interpolation of the values in table 4.7 of Ref. [TS], which is reproduced in table 
appendix |A} Then, we obtain 

^'wiT^'^wriw f " w>v " 

and S(Z, W) = 1 for W < Vq. As can be seen from figure [I] the effect of the screening correction itself is rather small 
and the theory behind it is well understood [T5], therefore we will not regard it a source of theoretical errors. 

Radiative corrections 

Radiative corrections are due to the emission of virtual and real photons by the charged particles present in (3- 
decay and it is crucial to correctly account for both the virtual and real photon contribution to ensure the proper 
cancellation of divergences. The finite part of the virtual photon exchange contribution does depend on the details 
of nuclear structure, but only affects the decay width and not the spectral shape. Therefore, the shape effect can 
be accurately computed. In our presentation we will suppress all terms which only affect the decay width, i.e. are 
independent of energy, and keep only the energy dependent parts. The radiative correction at order a has been 
computed by Sirlin [TH] for the /3-spectrum, with = p/W 

9, = 31nM„ - \ + 4 (^11) - U^2(W - W)] \ + %J ^ 



4 V P J \ 3W 2 n J \l + / 

+ ^ tanh- 1 (2(1 + (3 2 ) + ^""f^ " ^ann" 1 fl) , (12) 

where L(x) is the Spence function, defined as L(x) — dt/t ln(l — t). The complete correction is then given by 

G p {Z,W) = X + ^g fi . (13) 

The radiative correction for neutrinos, G„, has been derived in Ref. |20j and we apply the form given in Ref. |20j to 
the neutrino spectrum. The expressions are somewhat lengthy and involve integrals which can not be solved in closed 
form, therefore we do not quote them here. More recently, the radiative correction for neutrinos has been derived in 
closed form by Sirlin pTlJ^] and we find that the energy dependence of the new calculation, is numerically within 5% 
of the previous one [20]. Due to the fact that all spectra are normalized to unity, differences in terms which do not 
depend on energy will cancel. Thus, for our purposes the two results are equivalent and since the new result [21j is 
rather compact, we quote it here 

h v = 3\nM N + =f - [ I +8 ( 13 - 1 ) \n(2W$) 




tanh" 1 /? f 7 + 3f3 2 „ , \ 
+4 ^(_^-2tanh 1 \ , (14) 

where — p/E and p = \/w 2 — 1 and W — Wq — W v , where W v is the energy of the neutrino. The neutrino radiative 
correction then becomes, G V (Z,W V ) = 1 + a/(2-K)h v . Note, that the radiative correction in the /3-spectrum is quite 
large, whereas the correction for the neutrino spectrum is very much smaller. 



I would like to thank T. Schwetz-Mangold for bringing this new work to my attention. 
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Weak magnetism 



The term "weak magnetism" was coined by Gell-Mann and he proposed it as a sensitive test of the, then emerging, 
V-A theory of weak interactions 22J . The weak magnetism contribution to /3-decay belongs to the class of so called 
induced currents, i.e. currents which do not correspond to couplings present in the initial Hamiltonian and they 
appear only at finite momentum transfer. Out of these induced currents, the weak magnetism term typically yields 
the largest contribution to the shape of the /3-spectrum, see e.g. [23], and therefore we will only discuss this term. 
The effect of weak magnetism is to modify the spectral shape of /3-decay by a factor 

4 b 

1 + ^wmW with S WM = o TT~ m e > ( 15 ) 

— w 

where c is the Gamow- Teller matrix element, b — %/2/i, with fi the magnetic transition moment and M = AMn 
the mass of the nucleus. Experimentally, the transition magnetic moment or b can be determined by measuring the 
magnetic dipolc Ml 7 decay width, Tmi, of the corresponding isovector transition of the isobaric analog state and b 
is then given by 

This derivation relies on the concept of conserved vector currents (CVC), which is essentially a necessary result of 
gauge invariance. In the so called impulse approximation, one further assumes that the transition magnetic moment 
is entirely given by the intrinsic anomalous magnetic moments of the proton, //„, and neutron, fi n , and one obtains 
for w, the following simple expression 



W — ~^T7 

M N 



C 



v 



C A 



(17) 



Note, that in general there also would be a contribution to b (or fi) from the orbital angular momentum of the decaying 
nucleon, which we neglect here. Within the impulse approximation and neglecting orbital angular momentunj^j there 
is no dependence on nuclear structure and therefore, w has the same value for all Gamow- Teller /3-decays. This 
approximation is the basis for the statement in Ref. [23] that 4J| = |w = 0.5% MeV -1 is the universal value for the 
weak magnetism slope parameter in all /3 decays. The validity of the necessary approximation seems to be supported 
by the measured value for 4J| in the A = 12 iso-triplet of 12 B, 12 C and 12 N [2"5] . 

At the same time, one can use the CVC hypothesis and use equation[l6]to infer the weak magnetism slope parameter 
for a number of isotopes from their measured 7-decay parameters. In Ref. |26j this analysis has been performed on 
a set of 12 nuclei and a large variation of both 4^ and b was found, see table 11 of Ref. [33]. In table [i] we present 
a somewhat enlarged set of nuclei for which we could identify the necessary information from the Evaluated Nuclear 
Structure Data File (ENSDF) database [27]. The analysis follows the one outlined in Ref. [26]. For some of the nuclei 
we find slightly different results in comparison to Ref. [2B], which can be traced to changes in the nuclear database 
and we have, in all cases, provided the corresponding reference. 

This sample of nuclei is small, 13 nuclei only, and the majority is much lighter than typical fission fragments, which 
raises the issue of how representative this sample is for fission fragments. On the other hand, these nuclei summarize 
our current experimental understanding of weak magnetism. Setting these doubts aside, we can compute the average 
value of 4J| and its standard deviation as a measure of the error. In table [i] we have separated three isotopes based 
on their unusually large log ft values. In those cases, the Gamow- Teller matrix element is very small and thus the 
relative size of the weak magnetism contribution is enhanced. As a result, all of these 3 isotopes have large values 
of 4g. Obviously, these 3 large-log ft isotopes will completely dominate the average and standard deviation and 
therefore, here the question of how representative our sample is for fission fragments becomes acute. For the main 
analysis presented in this paper, we therefore will exclude these 3 isotopes and this yields the following mean and 
standard deviations 

^ = (0.67 ±0.26)% MeV" 1 , (18) 



5 Under these assumptions, there is also no contribution of the weak electric form factor, called d in the notation of 23 , which could be 
sizable. 
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8.2 71.8 
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0.646 


[28] 


12 B - 


> 12 C 


1+ 


->0+ 
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43.6 37.9 


11640. 
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0.62 




12 N - 
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1+ 
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43.6 37.9 
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m 


18 Ne - 
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m 
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0.046 129. 
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0.3 20.8 
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0.3 20.8 
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2.53 
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♦ 14 N 


0+ 
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0.0067 9.16 


1.096 x 10 9 


0.00237 


276. 


37.6 
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♦ 14 N 


0+ 




2313 


0.0067 9.16 


1.901 x 10 7 


0.018 


36.4 
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m 


32p _ 


> 32 S 


1+ 


-^0+ 


7002 


0.3 26.6 7.943 x 10 7 


0.00879 


94.4 


12.9 


mi 



TABLE I. Gamow- Teller decays and the associated parameters needed for a computation of the weak magnetism slope 
parameter using the CVC hypothesis. 



which is very close to the value obtained from using the impulse approximation, but we also see that there is significant 
variance. We will use the standard value of Ref. [23] and assign a 100% error to it. 
If we include the large-log ft isotopes, the corresponding result would be 

^ == (4.78 ±10.5)% MeV" 1 , (19) 
aE 

that is, a ten times larger mean with a many times larger relative error. A shift of 4-^ by +0.5 % MeV -1 causes a 



shift of the neutrino rate of about —1%. Thus, if one were to use equation 19 the resulting rate uncertainty would be 
about 20% and the neutrino flux shift found in Ref. [6j would be quite easily absorbed in this error bar. Even just a 
shift of the central value is sufficient to bring the neutrino fluxes nearly back to the original level. Obviously, large 
values of log ft imply relatively long half- lives and long-lived isotopes play only a small role in neutrino emission^] from 
fission fragments. However, if the endpoint energy is large, the relevant Gamow- Teller matrix element can be quite 
small without leading to a very long lifetime. Also, it is not clear how large weak magnetism would be in forbidden 
decays, since there the leading order nuclear matrix element is already quite small and thus all corrections tend to 
have relatively more weight. Thus, the size of weak magnetism, or more generally induced currents, is the one major 
source of theory uncertainty. In summary, the reactor anomaly |10j can be either attributed to new physics in the 
form of a sterile neutrino or to a significant extent to some not well understood nuclear physics. In particular, a 
detailed study of the breakdown of the impulse approximation in large-log ft nuclei in the relevant mass range would 
be quite helpful, but is beyond the scope of this paper. 



III. EXTRACTION OF THE ASSOCIATED NEUTRINO SPECTRUM 



For a single /3-decay branch the inversion from the /3-spectrum into the corresponding neutrino spectrum is straight- 
forward due to energy conservation, if one neglects recoil effects, which are O(E /(AM N )) ~ 10~ 4 . In particular, if 
the /3-spectrum has been measured with good precision, the neutrino spectrum can be known to the same level of 
precision without any reliance on theory. In practice, however, most of the nuclei involved have many /3-branches, 
most of which are identified by 7-ray spectroscopy and not by direct measurement of the /3-spectrum; in these cases, 
/3-decay theory as outlined in section |Tl] plays a crucial role in obtaining the associated neutrino spectrum. Of course, 
in a nuclear reactor a very large number of different nuclei and /3-branches contributes to the neutrino spectrum, 
which makes the task of a first principles calculation of the neutrino spectrum challenging and as we have describe 
in section [TTJ theory errors can play an important role, even if one were to assume that the actual nuclear data itself 
were flawless. The problems of this direct or a priori approach have been explored in detail in Ref. [B], with the 
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FIG. 1. (Color online) Shown is the relative size of the various corrections listed in equation [4] for a hypothetical /3-decay with 
Z — 46, A = 117 and Eo = 10 MeV. The upper panel shows the effect on the neutrino spectrum, whereas the lower panels 
shows the effect on the /3-spectrum. 

conclusion that the inevitable incompleteness of the nuclear databases makes it impossible to account for all neutrinos 
and reliance on measured total /3-spectra can not be reduced beyond a certain level and eventually virtual /3-branchcs 
have to be used in addition to the a priori computed spectra. 

Neutrino spectra from reactor neutrinos have been obtained by theory calculation by many authors [38 43 and for 
238 U this is still the only possibility^ The neutrino spectra, which were the state of the art prior to Ref. [6], were 
obtained by the exclusive use of virtual /3-branches, which parameters were determined by a fit to measured total 
/3-spectra [7HH]- This procedure is not independent of input from nuclear databases, either. Here, the information 
from nuclear databases enters in the form of the empirical mean proton number of the fission fragments as a function 
of Eo, which we call the effective nuclear charge, Z, and is computed according to 



Z(E ) 



f °° dE' V (E>) 



(20) 



with A being an appropriately chosen energy interval. 

From a mathematical point of view, the problem at hand is an inversion or unfolding problem: Given the measured 
total /3-spectrum Afp, one tries to infer the underlying distribution of /3-branches rj, which we take to be a continuous 
distribution over all possible endpoint^J Eq 

ATp(E e ) = J dE o N (E e ,E o ;Z)T)(E o ). (21) 

If we neglect the dependence on Z, problems like this one are known as Fredholm integral equation of the first kind. 
These problems are ill-posed and in general have no unique solution, therefore it is necessary to provide an additional 
constraint, a so called regularization scheme. The choice of regularization scheme is essentially an art and not a science, 
since one has to find the appropriate compromise between finding a reasonably stable solution and the introduction 



7 According to K. Schreckenbach a measurement of the 238 U /3-spectrum is currently undertaken at the FRM II in Garching. 

8 For the moment we neglect the issue of forbidden decays. 
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FIG. 2. (Color online) The basic inversion procedure, shown for a synthetic data set for 235 U. The green (thin, gray) line 
shows the neutrino residuals in 50 keV bins and the blue (thin, black) line shows the /3-residuals in 50 keV bins. The thick red 
line depicts the neutrino residuals in 250 keV bins. 

of an unacceptable bias into the result; for an introduction into the topic of regularized inversion see e.g. Ref. [44] . 
The regularization scheme, which has been developed for our case, relies essentially on averaging the solution over 
a finite range of E v , which corresponds to a top-hat filter. This procedure has been studied in great detail in [45] 
based on a synthetic data set. The use of a synthetic data set as proxy for real data allows to quantify the size of 
the bias and thus serves as a valuable tool to fine tune the regularization scheme. In Ref. [5J it was shown that a 
/3-spectrum computed from cumulative fission decay yields and the ENSDF database reproduces the measured total 
/3-spectrum for 235 U within an error of about 10% and therefore can serve as an accurate proxy for the actual data 
set. We will use a total /3-spectrum derived from the endpoints contained in the ENSDF databas^Jand we use fission 
yields for fission by 25meV, i.e. thermal, neutrons from the JEFF database, version 3.1.1 [H]. Note, that we neglect 
the contribution from epithermal and fast neutrons which would change the total spectrum by only a few percent 
and therefore, does not affect its suitability or accuracy as proxy. We also do not apply the small correction for finite 
irradiation time, this correction can be found in table VII of Ref. [BJ. In computing the resulting /3-spectra we apply 
the results obtained in section [n] We use this data set, which consists of approximately 550 isotopes and contains 
about 8 000 individual /3-branches, to perform a critical study of an inversion purely based on virtual branches. 

In Ref. [33] it was shown that given a precise enough measurement of the total /3-spectrum and accurate knowledge of 
Z the true neutrino spectrum could be recovered to better than within 1%. Unfortunately, several of the assumptions 
made in this analysis do not apply here: first, the available data was taken with a resolution of 50keV and secondly Z 
can be known only approximately, since the underlying database is incomplete, otherwise there would be no need for 
the use of virtual branches. Also, the analysis in Ref. [25] does not include the various corrections to a allowed /3-decay 
shape and the effect of statistical fluctuations of the actual /3-data was not taken into account. In the following we 
will quantify the impact of this various factors. 

A. Bias &: Statistical Errors 

The basic procedure is to take the measured /3-spectrum, for which, at this stage, we will substitute our simulated 
/3-spectrum and to start from the highest energy data point rig , where N is the number of data points. Then, one 

fixes a certain size s of a slice S and takes the data points Si = {n^~ s , n^~ s+1 , . . . ,n^}. We then fit an allowed 
/3-spectrum with a free endpoint E± and amplitude a\ to the data contained in S\. We continue this /3-branch to all 
energies and subtract it from the remaining N ~ s data points. This procedure is repeated till all data points have 
been fitted and as result we have a set of v = N/s endpoints and amplitudes {Ei, a^}, where v is the number of virtual 
branches. The set {Ei,ai} is a discrete approximation to T)(Eq) and can be used to compute the neutrino spectrum, 
by trivially inverting each virtual /3-branch into its corresponding neutrino spectrum. The result of this procedure for 



The endpoints, branching fractions, the assignment of the degree of forbiddeness and simple /3-spectra which were computed according 
to Ref. j6], i.e. without the terms C, S, G„, were kindly provided in machine readable format by the authors of Ref. [6]. 
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FIG. 3. (Color online) The neutrino flux from the inversion of 1000 random realizations of a synthetic /3-spectrum for 235 U 
relative to the mean outcome of these 1000 trials. The red (thin, gray) lines show a subset of 100 trials. The thick, ragged 
black line shows one particular example. The smooth thick, blue (black) lines show the standard deviation in each bin, which 
is the same as the square root of the diagonal elements of the covariance matrix. The dark green (thick, dark gray) lines are 
the statistical error of the /3-spectrum scaled up by a factor of 7. 



235 U with "data" in bins of 50 keV and s = 5, corresponding to 30 virtual branches, is shown as green (thin, gray) 
line in figure [2j The thin blue (black) line shows the /3-residuals, which are obviously very small. 

The 50keV anti-neutrino residuals show oscillations which are due to the simple fact that the /3-spectrum is non- 
zero at the endpoint. The amplitude of these oscillations in increased by the inversion procedure and this increase 
is a direct consequence of the problem being ill-posed in the mathematical sense. To obtain a more regular solution, 
we can average our result and obtain the thick red line by averaging over 250 keV bins. Clearly, this regularized 
solution shows oscillations of much smaller amplitude, especially in the region of interest between 2 and 8 MeV. We 
also observe that the remaining oscillations do not average to zero, e.g. between 4 and 6 MeV both, the regularized 
and non-regularized, solutions predict a neutrino flux which is too low by about 1%; the result is biased. A finite bias 
is not surprising and quite common in this type of inversion problem. Using a synthetic data set it is straightforward 
to quantify the bias by comparing the result of the inversion with the known true flux, as done in figure [2] We, also, 
can vary the starting point of the inversion by removing an increasing number of consecutive data data points at the 
high energy end and repeat the procedure. The overall features stay the same, but the phase of the oscillations of 
the non-regularized solution changes gradually and once we have removed s data points we obtain the initial solution 
again. This variation of the phase introduces a slight spread in the magnitude of the bias, which we will assign as 
error of the bias correction. 

The oscillations apparent in figure [2] show that the inversion procedure amplifies small perturbations and may be 
not stable with respect to statistical fluctuations, which in real data are inevitable. Given the statistical errors of the 
/3-spectra we can add random fluctuations with the same variance to our synthetic /3-spectra and for each realization 
of a randomized /3-spectrum we can perform an inversion. We have done this for 1000 random data sets and the 
results are shown as red (thin, gray) lines in figure [3j where we plot the relative deviation from the mean. The thick 
black lines shows one characteristic example. It is quite obvious that there are very strong anti-correlations between 
neighboring bins. These anti-correlations are due to oscillations in the inversion procedure, which are excited by 
statistical fluctuations. Since, these oscillations are unphysical and merely a feature of the inversion procedure, we 
do not want to propagate these oscillations into the final result. The thick smooth, black line shows the standard 
deviation in each bin, which turns out to be the same as the square root of the corresponding diagonal element of the 
covariance matrix of the 1000 randomized inversions. The range delimited by the standard deviation is a good estimate 
of the statistical error, therefore, we will drop all off-diagonal elements of the covariance matrix. The resulting errors 
are found to be several times larger in the inverted neutrino spectrum than they are in the underlying /3-spectrum. 
The errors of the /3-spectrum are shown as dark green (thick, gray) lines but multiplied by a factor of 7. Note, that 
this scaling factor is not universal and depends sensitively on the details of the inversion, like bin width. 
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FIG. 4. (Color online) The effective nuclear charge Z of the fission fragments of 235 U as a function of Eq. The area of the 
each box is proportional to the contribution of that particular Z to the fission yield in that energy bin. The lines are fits of 
quadratic polynomials: black - ENSDF database, blue (dark gray) - adding those isotopes missing from ENSDF by assuming 
that there is only one /3-branch, each with Eo = Qp, red (light gray) - maximum pandemonium as defined in the text. The 
blue (dark gray, rightmost), filled boxes show the resulting distribution of adding the missing isotopes with Eq = Qp. The red 
(light gray, leftmost), empty boxes show the distribution in the maximum pandemonium approximation. 



B. Effective average nuclear charge 

Finally, we study the degree to which Z can be determined and how the incompleteness of nuclear databases affects 
the inversion procedure. We fit a second order polynomial in Eq, where we weight each /3-decay branch by it fission 
yield Ya,z and the corresponding branching ratio b\ z , the result is shown as a black line in figure 4 The black boxes 
are a two-dimensional histogram of the Z distribution of the fission fragments using the ENSDF database. Note, 
that this distribution is bi-modal, because fission proceeds into two asymmetric fragments, and it is therefore, quite 
unexpected that a simple energy-dependent mean (Z) actually is sufficient to compute accurate neutrino spectra, as 
previously shown in [45) and in figure [2j 

We know from comparing the fission yields with the ENSDF database how many and which nuclei are missing 
and although we have no detailed information on /3-branches for those nuclei, ENSDF does contain the Qp for all of 
the missing nuclei. Thus, we can try to bracket the effect of the missing nuclei by computing Z under two different 
approximations. First, we assume there is indeed only one /3-branch with Eq = Qp for each missing isotope. In this 
case, the effect of the missing isotopes is to add to the high energy part and increase Z by a rather small amount, as 
shown by the blue (dark gray) curve in figure |4j In nuclei far from stability, there will be very many, but individually 
quite weak /3-branches, which makes it very difficult to infer the correct level scheme and /3-decay branching fractions 
from 7-spectroscopy, this is termed as pandemonium effect |47| . In Ref. [5] the pandemonium effect is taken care 
of by replacing about 50-100 of the nuclei in the ENSDF database by data which was obtained specifically avoiding 
the pitfalls of the pandemonium. Here, we would like to use an approximation to this method, which we will call 
the "maximum pandemonium". We take the 17.5% most neutron rich nuclej^] from the ENSDF database and all 
the Qp values of the missing nuclei and replace all their /3-spectra using the following algorithm: We distribute the 
/3-decay strength evenly between and Qp using about 10 decay branches. The result is shown as red (light gray) 
line in figure |4j This should yield a reasonable approximation, at least when one averages over a large number of 
nuclei, as we presently do. The maximum pandemonium approximation is of course rather crude and would be not 
sufficient for a direct a priori calculation of the neutrino spectrum, but it allows us to gauge the magnitude of the 
effect of the incompleteness of the nuclear database we use. We repeat this exercise for the other isotopes as well and 
the results on Z are summarized in table |ll| The induced systematic error is correlated between the isotopes since it 
has a common physical cause. Note, that this is the only place at which information from nuclear databases enters 
the computation of the neutrino spectrum using only virtual /3-branches. In principle, the nuclear mass A will show 



17.5% was chosen because the resulting deficit in the total /3-spectra is very similar in magnitude and shape to the one shown in figure 5 
of Ref. [§]. 
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Isotope 


CO 


Cl 


C2 


235jj 
239p u 

241 Pu 


48.992±°, 164 
49.650l°. 214 
49.906l°. 178 


-0.399l° 161 
_ 447 +o.o36 

-0.510±° 160 


-0.084l° 044 
-0.089l°. 016 
-0.044l2. 052 



TABLE II. The coefficients for the parametrization of Z for the various isotopes. Also given are the range as obtained from 
the maximum pandemonium approximation. Ci is the coefficient of Eq in a second order polynomial in Eq. Note, that in our 
definition Z is the charge of the parent nucleus. 

a similar change as a function of Eq, but variations of A have a much smaller impact on the neutrino spectrum and 
therefore fixing the effective nuclear mass, A, at A = 117 does not lead to any relevant additional bias or error. The 
value for A is chosen based on the average nuclear radius, where the average is determined by weighting each isotope 
with its cumulative fission yield and we find a range of about 116-118 for 235 U to 241 Pu. 



IV. APPLICATION TO 235 U, 239 PU AND 241 PU 



In this section, we will apply the results obtained so far to a direct inversion of the /3-spectra which were measured 
in the 1980s at ILL. We use the /3-spectra for 235 U as presented in Ref. [7,, for 239 Pu from Ref. [HHf^and for 241 Pu 
from Ref. [9]. The original data for all 3 isotopes was recorded in 50keV bins, but unfortunately was published only 
in 250 keV bins, which is not sufficient for our purposes. We, however, were able to obtain finer binned data from 
one of the authors [IS] of Refs. [7HS]: for 235 U we have data in 50keV bins from 1.5 to 9.5 MeV, for 239 Pu we have 
data in 100 keV bins from 1 to 8 MeV and for 241 Pu we have data in 100 keV bins from 1.5 to 9 MeV. Thus, we have 
data for all 3 isotopes only in the energy interval between 1.5 and 8 MeV and, since the threshold for inverse /3-decay 
is 1.8 MeV and the results in Ref. [B] range from 2 to 8 MeV, we will present result only for this interval. Together 
with the original data, also came the size of the statistical error in each bin, which when summed up over 250 keV 
bins agree with the values published in Refs. [7-9], as do the data itself. Note, that the data on 235 U not only has the 
smallest bin size but also the smallest statistical errors, whereas 239 Pu has the largest errors. There is a calibration 
or normalization error for each data set and we take these values directly from Refs. [7H3] and list them in column 9 
of tables IVIlllXl 



For each isotope we perform the inversion as outlined in section III but now for the actual data, we use 30 virtual 
branches with a slice size of 5 for 235 U and 23 virtual branches for 23y Pu and 25 for 241 Pu both with a slice size of 3. 
We have tested that variations of the slice size around the values given only mildly affect the results and the chosen 
values provide the smallest errors around 4 MeV. Also, with these choices the energy range of a slice is about the 
same in all 3 isotopes and is close to 250 keV. By using synthetic data sets for each isotope we determine the bias 
and statistical errors, the later is shown column 5 of tables |VII|IX| where the square root of the diagonal elements 
of the covariance matrix are given. We apply the resulting bias correction, column 3 of tables |VIipX| to the results 
extracted from the actual data and assign the spread of biases obtained by varying the starting point of the inversions 
as error, which is given in column 6 of tables [VlT|lX| s ee als o section [ill A| Then we perform an inversion using the 



upper and lower ends of the ranges for Z from section IIIB and for the weak magnetism correction from section |l"l| 



The relative differences to the initial inversion result are listed as errors in columns 7 and 8 of tables |VIipX| Finally 
the error components in columns 5-9 are added in quadrature and quoted as total error in column 10 of tables [Vn|lX| 

The total error quoted in tables VII||lX is only an indication of the actual errors since there are correlations between 
the various bins. The statistical error is uncorrelated between bins as well as the bias error. Both the errors due to 
Z and weak magnetism (WM) are fully correlated between all bins and all isotopes. The normalization errors is fully 
correlated between bins and between isotopes, since they same apparatus and method to determine the normalization 
was used for all 3 isotopes. Due to the asymmetric nature of the errors from Z and weak magnetism (WM) they can 
not be included exactly in a total covariance matrix. However, already the combined errors form Z and WM are close 
to symmetric and the overall errors are very close to symmetric, thus it is a reasonable approximation to neglect the 
small asymmetry in errors for most applications. 

For completeness, we also provide a parametrization of our results using the familiar exponential of a 5 th order 



11 The actual data on electron spectra is the same in both sources, however the neutrino spectrum was re-derived in the later reference 
using updated nuclear databases. 
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FIG. 5. (Color online) Comparison of our result for 235 TJ with previous inversions, labeled ILL for the results from Ref. [7] 
and labeled 1101.2663 for the results from Ref. [6]. The thin error bars show the theory errors from the effective nuclear charge 
Z and weak magnetism. The thick error bars are the statistical errors, whereas the light gray boxes show the error from the 
applied bias correction. The green line, referred to as simple, shows the result, if we use the same description of /3-decay as in 
Ref. j6]. The black line, referred to as ILL inversion, shows our result if we completely follow the procedure outlined in Ref. [7], 
including their effective nuclear charge. 



polynomial [ol 135]: 



b v {E v ) =exp l^aiBt 1 



(22) 



In performing the fit we follow the description given in Refs. [5] 3PJ but we do not include contributions from the 
error on Z and WM since they are correlated between isotopes a nd d o not change the fit appreciably. The resulting 
best fit parameters and the minimum \ 2 values are given in table III 



Isotope 


Xmin 


Ctl Q2 C*3 Ct4 &5 C(6 


235jj 


49.3 


4.367 -4.577 2.100 -5.294 x 10" 1 6.186 x 10" 2 -2.777 x 10 -3 


239p u 


20.8 


4.757 -5.392 2.563 -6.596 x 10" 1 7.820 x 10" 2 -3.536 x 10" 3 


241 Pu 


15. 


2.990 -2.882 1.278 -3.343 x 10" 1 3.905 x 10 -2 -1.754 x 10 -3 



TABLE III. Result of a fit of a 5* order polynomial to the logarithm of the flux. The number of degrees of freedom is 25 — 6. 

Obviously, the fit for 235 U is quite bad, with a x 2 /dof of more than 2. Also, for all 3 isotopes the fit parameters 
are highly correlated and we therefore do not provide any fit errors or correlation matrices, since it seems doubtful 
whether these could be used to model the errors in the underlying neutrino fluxes. This parametrization should not 
be used for actual data analysis or error propagatiorp^] Instead, we recommend to rebin our results by using linear 
interpolation and integrating the resulting fluxes over the new bins. In the same fashion the errors can be rebinned 
and we provide results for 250 keV and 50 keV bins in machine readable format [5D] . 

V. DISCUSSION 

The /3-spectra from Refs. [THPJ have been previously inverted into neutrino spectra and we, therefore, will start by 
comparing our result with previous results [HHPJ- In figure [5] we show our result (thick blue/black line) relative to the 
results presented in Ref. [7j, which is denoted as We clearly observe that our results point to a significantly 



12 Nonetheless, this parametrization can be safely used to extrapolate our results to higher energies, since there, the errors of the actual 
fluxes are sufficiently large to render the inaccuracies of this parametrization harmless. 
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enhanced flux with respect to </>ill- The results from Ref. [BJ are shown as red (thin, dark gray) line and up to about 
5.5 MeV there is good agreement, whereas at higher energies the result of Ref. [6J exhibits strong oscillation relative to 
both </> IL l and our result, which may be due the coarser binning of the /3-spectra used in Ref. [BJ. The thin error bars 
are the sum of the theory errors, i.e. the effective nuclear charge Z and weak magnetism WM, the thick error bars 
show the statistical errors and finally the gray boxes show the error from the bias correction. The theory errors are 
correlated between all bins, the statistical and bias correction error are uncorrelated. Not shown is the, common to all 
data sets, normalization error. The thin green line shows our result, if we use the same, simplified, description of the 
/3-decay spectrum as in Ref. jS] and we see that it is responsible for a sizable fraction of the difference of our result with 
the others at energies above 5 MeV. Finally, the thin black line shows our result, if we follow the procedure outlined 
in Ref. [7], which in particular includes their form of Z and the fact, that all corrections to a allowed /3-spectrum 
shape are applied in an average fashion and we recover the original result quite well. Similar figures can be made 
for the other isotopes and the results, generally, would be similar, too; the error bars would be much larger due to 
the 100 keV binning of the original data and the larger statistical errors in the /3-spectra. As far as the origin of the 
upward shift in flux relative to the original results [7H9] is concerned, we agree at some level with the comments made 
in Ref. [BJ , where at energies below 4 MeV the main effect is ascribed to the branch-by-branch implementation of the 
various corrections to a allowed /3-spectrum shape. We would like to point out that both an implementation in each 
virtual branch as well as an implementation in each physical branch in an a priori calculation yield the same effect. 
Moreover, it is likely that a newly derived average correction based on modern nuclear databases maybe nearly as 
accurate as a branch-by-branch implementation. Of course, to derive this new average correction would presumably 
imply to actually study the corrections first branch-by-branch, so nothing would be gained. At energies above 5 MeV, 
the new form of the effective nuclear charge Z as given in table [n] does make up a large fraction of the shift. Note, that 
in contrast to Ref. [BJ we do find that an effective nuclear charge is sufficiently accurate and the effect of the actual 
distribution of Z can be approximated quite well by its, properly defined, average; this conclusion is supported by the 
results of Ref. [45 . At high energies the contribution from the wave function convolution and screening correction on 
the /3-shape contribute significantly to the overall shift. Once we combine all factors we find about the same upward 
shift in all isotopes, however we would not characterize it as constant over all energies. 

In table [IV] we quantify the mutual agreement or disagreement of the various fluxes and associated inverse /3-decay 
event rates, where we use the cross section from [5T]. Inverse /3-decay is frequently used for the detection of reactor 
anti-neutrinos. In this table we list for each isotope the total difference over all energy bins from 2 — 8 MeV of the 
fluxes from Refs. [7HS] labeled 0ill an d the ones from Ref. [BJ labeled </>mlf to our results. We provide the total error 
as well as the individual contributions from statistics, theory and the bias correction. The ratio of the total difference 
to the total error is what we call the "rate significance", i.e. the significance of the total flux change. We find the 
total flux change relative to c/>tll to be in the range +(2.4 — 3.2)%. The upward shift in fluxes is significant at the 
(3.7 — 4.0) a level. This shift is slightly larger than the one found in Ref. [BJ, however if we compare our total flux 
with respect to 0mlf we find very small differences for 235 U and 241 Pu. For 239 Pu the situation is different and we 
have a lower total flux than </>mlf 7 which however is entirely due to the first energy bin and could be due a spurious 
oscillation in the inversion procedure in either analysis. If we perform the same comparisons for inverse /3-decay event 
rates, rows labeled as .Rill an d -RmlFj we find that in total event rates the upward shift with respect to -Rill is a 
bit larger than for the fluxes with +(3.7 — 4.7)%, which is due to the relatively large shift at high energies which is 
enhanced by the energy dependence of the cross section. These rates shifts are significant at (2.4 — 3.0) a and the 
reason for the significance being smaller than for the fluxes is again the larger weight of high energies which do have a 
larger shift but also larger errors. As far as total inverse /3-decay rates are concerned, the agreement with the results 
of Ref. [BJ is excellent and no difference reaches more than 1 a. 

In the last column of table IV we provide the "shape significance" , i.e. the square root of the x 2 -difference between 
the various flux and inverse /3-decay rates. For this calculation we assume a free normalization, which is equivalent to 
subtracting any rate differences. For the shape difference of fluxes we assume infinite statistics, whereas for the inverse 
/3-decay shape difference we assume 10 6 events. We see that shape differences can be large and quite significant at 
even more than 5 a. This is, however, based only on the flux errors and does not include further complications like 
experimental errors or errors on the isotope composition of the reactor and burn-up. Thus, at this level of analysis it 
remains unclear whether an actual neutrino measurement could distinguish between the various flux models. 

Apart from the quantitative difference we discussed previously, there are important qualitative differences of our 
result with respect to earlier inversions. We present, for the first time, a careful and detailed error analysis of both 
the underlying theory of /3-decays as well as of the inversion procedure itself. We found that most corrections to 
the statistical /3-decay shape, like finite size effects, are well understood and the associated errors are only a fraction 
of the correction itself, see section [TTJ However, we also found that the treatment of induced currents in /3-decay is 
subject to considerable theoretical uncertainty We studied in detail the example of weak magnetism as the largest 
contributor and found by comparison with data, the impulse approximation, which is the basis for assuming that the 
weak magnetism correction is universal for all beta decays, is apparently only accurate within about 100% for most 
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235-q 


shift 

% 


1 a errors [%] 
total stat theo. bias 


a significance 
rate shape 


(frlLL 

Rill 


2.4 
3.7 


0.7 0.1 0.6 0.1 
1.5 0.0 1.5 0.1 


3.7 4.5 
2.4 3.7 


0MLF 

Rmlf 


0.4 
1.2 


0.7 0.1 0.6 0.1 
1.5 0.0 1.5 0.1 


0.6 6.0 
0.8 5.4 


239p u 


shift 

% 


1 a errors [%] 
total stat theo. bias 


a significance 
rate shape 


4>ILL 

Rill 


2.9 
4.2 


0.7 0.3 0.6 0.2 
1.5 0.4 1.4 0.2 


3.9 3.2 
2.8 3.2 


0MLF 

-Rmlf 


^ 

-°tl 
1.4 


0.7 0.3 -0.6 -0.2 
1.5 0.4 1.4 0.2 


0.5 4.7 
1.0 4.6 


241 Pu 


shift 

% 


1 a errors [%] 
total stat theo. bias 


a significance 
rate shape 


(frlLL 

Rill 


3.2 
4.7 


0.8 0.3 0.7 0.2 
1.6 0.2 1.5 0.2 


4.0 2.6 
3.0 2.6 


0MLF 

Rmlf 


0.4 
1.0 


0.8 0.3 0.7 0.2 
1.6 0.2 1.5 0.2 


0.5 3.7 
0.7 3.7 



a This difference and its sign are driven by the first bin from 1.875 — 2.125 MeV and thus the inverse /3-decay rates do not reflect this. 

TABLE IV. These tables compare our result in relation to the fluxes from Refs. 7-9 labeled as 4>ill or the inverse /3-decay 
event rates -Rill- The fluxes and event rates from Ref. [5] are referred to as </>mlf and -Rmlf- For the column labeled rate 
we do not include any normalization errors since they are common to the underlying data, for the column labeled shape we 
assume a free normalization. The significance in the shape column for event rates assumes 10 6 inverse /3-decay events. 

nuclei and fails completely for a small subset of large-log ft nuclei. For the analysis presented so far we assumed a 
100% error on the weak magnetism corrections, i.e. we neglect the large-log ft nuclei. If we assume that our limited 
sample of nuclei were representative for fission fragments this error would be 10 times larger and the sole dominating 
contribution to the overall error budget, in which case all flux models would be equivalent within these much larger 
error bars. Even if we assume further, that the error bars would remain small, merely adjusting the central value 
within this larger range to about 4 times its current value would allow to largely eliminate the upward shift and to 
bring both the calculation presented in Ref. jS] and our result into agreement with the previous ILL fluxes. As a 
consequence the new found support for sterile neutrinos as presented in Refs. |10|lll] would be weakened considerably. 

We also investigated in detail how errors propagate through the inversion procedure and quantified the bias in this 
method by using synthetic data sets. It is important to realize that the underlying mathematical problem is ill-posed 
and thus we are not dealing with a simple (or complex) fit of a model to data. Therefore, the use of synthetic data 
sets as first shown in Ref. [55] is essential to gain insight into error propagation and errors inherent to the method. 
The assumption that the statistical errors of the /3-spectrum would find a one-to-one correspondence in the resulting 
neutrino spectra is not supported by our Monte Carlo simulations. Also, the size of the bias is non-negligible and 
would be much larger if we had to use /3-spectra binned in 250 keV bins. The combination of bias and statistical 
inversion errors from our synthetic data sets yields the total conversion error. We also would like to emphasize that 
the /3-residuals obtained here are significantly smaller than in Ref. [5] and exceed nowhere the errors of the final result. 

In all methods of inversion, input from nuclear databases is required, since the measured /3-spectrum can be fit by 
a wide range of different branches by adjusting the nuclear charge. The problem is, that nuclear databases tend to be 
incomplete and may contain significant errors in level assignments, e.g. the pandemonium effect. There are different 
strategies to deal with these shortcomings and in Ref. [B] the attempt was made to address these shortcomings directly 
at the database level by collecting the best available experimental data and supplement missing data by appropriate 
theoretical models. This resulted in a database which is about 90% complete and the remainder was treated by 
inversion, i.e. with virtual branches. This method, if carefully applied, may have the potential to ultimately yield 
the most accurate results. Here, wc followed a different strategy which tries to minimize the use of nuclear databases 
as much as possible. In our inversion procedure which is solely based on virtual branches, the only place at which 
information from nuclear databases enters directly is in the form of the effective nuclear charge Z. Since Z is an 
average quantity, it is relatively straightforward to estimate the impact of both incompleteness, which we know from 
comparing fission yields with the known entries in the database, and of errors in data itself like the pandemonium 
effect. It turns out, that Z is amazingly robust with respect to these issues and never exceeds the error contributions 
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from nuclear theory. Therefore, we reason that even improved databases would barely change our result. 

In combination all our, quantified, error sources lead to total errors which are close to the ones found previously [BJ. 
At high energies our errors are consistently larger, whereas at lower energies they are either slightly smaller for the 
case of 235 U or slightly larger for the 2 Plutonium isotopes, where "slightly" means of the order 0.5%. 



VI. SUMMARY AND CONCLUSION 



We present a new inversion of the total /3-spectra from fission fragments obtained at ILL [7H5] employing only 
virtual /3-branches. Our analysis is based on the original data [Hj, which have finer binning, 50keV for 235 U and 
100 keV for 239 Pu and 241 Pu, compared to the published numbers [ZHS], but are otherwise identical. In comparison 
to Refs. [7HS] and [BJ we treat the theoretical /3-spectrum shape to a higher level of accuracy by including several 
sub-leading corrections, not previously considered in this context. Specifically, we include a more accurate treatment 
of finite size effects, screening and the radiative QED correction for the neutrino, for details see section [II] All of these 
corrections to a allowed /3-spectrum are applied at a branch- by-branch level, which was not done in Refs. [IrlS], where 
only an average correction was applied. The only point at which nuclear st ructur e information enters this calculation 



is in the form of the effective nuclear charge Z, which is derived in section IIIB Combining all these ingredients we 



find the fluxes as listed in tables VII||TX Our results averaged over all energies are in very good agreement with the 



results found in Ref. [6] and confirm the overall 2-3% shift relative to the original inversion [7H9], while at the same 
time the /3-spectrum residuals are very small, which was not the case in previous analyzes [BJ. 

We put a particular emphasis on evaluating the errors of the final neutrino fluxes. We analyzed the associated 
theory errors and found, that only the contributions from induced currents, for which we studied the weak magnetism 
term as the leading contribution, have appreciable theory errors. Our overall theory error is similar to the values 
previously quoted [6-9 , however the origin is different, in our case, it being solely due to the difficulties to estimate 
the size of the weak magnetism correction. We also point out, that there are certain nuclei like 14 C in which the 
weak magnetism correction is anomalously large due to a suppression of the Gamow- Teller matrix element. A small 
abundance of nuclei of this kind in fission fragments could enhance the weak magnetism term greatly, see section [IT] 
As a matter of fact, this enhancement could be large enough to account for a large fraction of the flux shift found here 
and in Ref. [BJ and thus may provide a Standard Model explanation of the reactor anomaly [TU] without the need for 
sterile neutrinos. While we are not advocating this solution to the reactor anomaly throughout this work, it seems 
that Occam's razor warrants a closer look into this possibility. We provide an estimate of the effects incompleteness 
and inevitable inaccuracies of the nuclear structure data have on the determination of the effective nuclear charge Z, 
see section |IIIB| We evaluate the inversion errors by using synthetic data sets based on cumulative fission yields and 
the ENSDF database and use Monte Carlo simulations to "measure" the bias and statistical error. We find, that the 
bias generated by our method is quite small, in agreement with earlier results [45] . The statistical errors are sizable 
and typically several times larger in the neutrino spectrum than in the underlying /3-spectrum, see section [HI A| The 
resulting total errors are similar to previous ones [BJ at low energies, however they are larger for the few highest 
energy bins. Based on these errors we find averaged over all energies excellent agreement with the results of Ref. [BJ, 
however if we look at the energy dependence we do find significant shape differences, particular at high energies. 
These difference should be accessible to measurement, provided a large enough neutrino event sample of the order 10 6 
events can be obtained, see table |IV| 



As far as the origin of the shift relative to the original results [THS] is concerned, we find that both the form of the 
effective nuclear charge Z and the use of an average correction to the allowed /3-shape contribute at a similar level. 
There are now two independent and complementary analyzes, which find the same sign and approximate magnitude 
of this shift. However, relevant differences between the two analyzes do exist, which can be briefly summarized as: 
our analysis uses finer binned input data and consequently has much smaller /3-spectrum residuals; we use a more 
detailed treatment of corrections to the allowed /3-spectrum shape and derive the associated errors from this more 
detailed treatment; we compute the inversion errors, both bias and statistical, based on synthetic data sets; finally, 
we estimate the errors due to incomplete or incorrect nuclear structure data. 

The aforementioned differences make it seem worthwhile to analyze existing reactor data using our results in a spirit 
similar to Ref. [10] and to investigate whether an actual neutrino spectrum measurement can distinguish between the 
two calculations [52]. The differences tend to be largest at high neutrino energies, where typically also the largest 
differences between the neutrino fluxes from the various isotopes are found [49] , thus the impact of these new fluxes 
and associated errors on neutrino safeguards schemes warrants further study. 
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Appendix A: Supporting Data 
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-11.223 
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23.9774 


-12.6534 
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«3 
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an 
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-471.299 


662.191 


-305.68 


as 


-1.5632 


33.4192 


-255.133 


938.53 


-1641.28 


1095.36 



TABLE V. A reproduction of table 1 of Ref. [T5] . 



Z 
N(Z) 



1 8 13 16 23 27 29 49 84 92 
1.000 1.420 1.484 1.497 1.52 1.544 1.561 1.637 1.838 1.907 



TABLE VI. A reproduction of table 4.7 of Ref. [16] 
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Appendix B: The Fluxes 
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TABLE VII. Results for the U anti-neutrino spectrum. This spectrum corresponds to 12 h irradiation time. The /3-residuals 
are given for information only and they do not enter into the computation of the errors. The errors in columns 5 and 6 are fully 
uncorrelated, whereas the errors in columns 7 and 8 are fully correlated between bins and isotopes. The errors in column 9 are 
fully correlated between bins and isotopes. 
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TABLE VIII. Results for the 239 Pu anti-neutrino spectrum. This spectrum corresponds to 36 h irradiation time. The /?- 
residuals are given for information only and they do not enter into the computation of the errors. The errors in columns 5 
and 6 are fully uncorrelated, whereas the errors in columns 7 and 8 are fully correlated between bins and isotopes. The errors 
in column 9 are fully correlated between bins and isotopes. 
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+0 

-3. 


+3. 

-2.6 


2. 


+9.1 

-9.4 


8 


0.7 


-3.5 


1.06 x 10" 3 


12. 


0.63 


+ 
-3.3 


+3.2 

-2.7 


2. 


+ 12. 
-13. 



TABLE IX. Results for the 241 Pu anti-neutrino spectrum. This spectrum corresponds to 43 h irradiation time. The /3-residuals 
are given for information only and they do not enter into the computation of the errors. The errors in columns 5 and 6 are fully 
uncorrelated, whereas the errors in columns 7 and 8 are fully correlated between bins and isotopes. The errors in column 9 arc 
fully correlated between bins and isotopes. 
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